Superconducting and pseudogap phases from scaling near a Van Hove singularity 
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We study the quantum corrections to the Fermi energy of a two-dimensional electron system, 
showing that it is attracted towards the Van Hove singularity for a certain range of doping levels. 
The scaling of the Fermi level allows to cure the infrared singularities left in the BCS channel 
after renormalization of the leading logarithm near the divergent density of states. A phase of d- 
wave superconductivity arises beyond the point of optimal doping corresponding to the peak of the 
superconducting instability. For lower doping levels, the condensation of particle-hole pairs due to 
the nesting of the saddle points takes over, leading to the opening of a gap for quasiparticles in the 
neighborhood of the singular points. 
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The study of the electronic properties of the cuprates 
represents nowadays a great challenge from the theoreti- 
cal point of view, as the phenomenology of these materi- 
als has become increasingly rich during the last decade. 
There have only been a few attempts to develop a theory 
that may encompass the main experimental features |Q , 
including the antiferromagnetism of the undoped com- 
pounds and the pseudogap phase above the supercon- 
ducting transition. Additionally other proposals have 
focused on the mechanism of superconductivity, stress- 
ing the role played by antiferromagnetic fluctuations ||] 
or by the proximity to a Van Hove singularity (VHS) in 
the doped materials ||. 

The later approach has received much attention re- 
cently, since it establishes a natural competition be- 
tween magnetic and superconducting instabilities in a 
two-dimensional (2D) system The investigation 

of the model of electrons near a VHS is delicate due to 
the appearance of logarithmic divergences in perturba- 
tion theory. In a renormalization group (RG) framework, 
one has to handle infrared singularities which arise after 
renormalizing away the leading logarithm, as the energy 
dependence of some quantities comes in powers of a log- 
arithm square Q . 

Most part of the analyses of the problem have been 
made fixing the Fermi level at the VHS from the start. 
This questions the naturalness of the predicted instabil- 
ities, that rely critically on the proximity to the singular 
density of states. The Fermi energy is actually a dynam- 
ical quantity that is shifted by quantum corrections. It 
has been shown that the VHS has the tendency to at- 
tract the Fermi level of the electron system |T^-[l3| . It is 
therefore more appropriate to let the chemical potential 
free to evolve as the states are integrated in the quantum 
theory. This also solves at once the problem of the in- 
frared divergences, as the shift of the chemical potential 
from the VHS acts as an infrared cutoff in the logarithmic 
dependences left in the renormalization. 

We illustrate the above ideas in the case of the t — 
t' Hubbard model, which has a dispersion relation with 
saddle points at A = (tt, 0) and B = (0, it). We deal with 
a wilsonian RG approach in which the chemical potential 



/i is originally placed away from the VHS, and electron 
modes in two thin slices about energies fi + A and fi — A 
are progressively integrated out |Q. At each RG step, 
the chemical potential is free to reaccomodate due to the 
self-energy corrections from the charge integrated out. 

The behavior of \x as A — > can be obtained by solv- 
ing the frequency and momentum-independent part of 
the Schwinger-Dyson equation 1/G = 1/Gq — E. The 
charge in the slice dA is given in terms of the density 
of states n(e) by n(/i — A)rfA. It couples through the 
forward-scattering vertex F in the usual Hartree and 
exchange diagrams. To improve the perturbative ap- 
proach, one has to take into account that F is also a 
scale-dependent quantity whose value vanishes in 
the proximity to the VHS according to the expression 
F{e) w F /(l - J F 1 log(e)/(47r 2 t)). Thus, upon integra- 
tion of high-energy modes in a thin slice dA, the chemical 
potential /i is shifted according to 



^ = F( M /( M -A))n(/i-A) 



(1) 



The differential renormalization given by Eq. (|l|) leads 
to the scaling of /j, as a function of A. It can be seen that, 
in a certain range, the flow of the chemical potential is 
attracted by the VHS as A — > 0. If F were constant, the 
fixed-point condition —1 + Fn(p) = 0, with p, = /j, — A, 
would give /i(A). Taking into account the scaling of the 
F vertex, it is found that the renormalized value of /x lies 
very close to to the level of the singularity, for appropriate 
values of the bare chemical potential. 

We have represented in Fig. |]the results of solving the 
scaling equation with a model density of states n(e) — 
\og(t j e) / '(4ir 2 t) for e < 0.5t, and constant elsewhere. 
We observe that there is a range of nominal values of the 
chemical potential in which this is attracted towards the 
VHS. As a consequence of that, there is a range of filling 
levels that are forbidden above the level of the singularity. 

In an open system, the pinning of the Fermi level to 
the VHS can be realized by exchanging particles with the 
environment, as greater stability is attained then in the 
system |13| . In a closed system, however, the mechanism 
of pinning can be only realized if the Fermi line deforms 
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to become closer to the saddle points of the dispersion 
relation. This is consistent with the renormalization of 
the Fermi line observed in numerical studies of the Hub- 
bard model, which is translated into a shift of the next- 
to-nearest- neighbor hopping t' Such an effect is 
also in agreement with recent experimental observations 
reported in Ref. p7[ . It has been shown there that the 
Fermi line of Bi-2212 bends progressively as the doping 
level is increased, so that even in the overdoped regime 
it does not lose its hole-like character. 




vant contributions ~ (dA) 2 at each RG step 

The converse statement is also true, and the vertices 
with the BCS kinematics Vi n t ra and V um k only get renor- 
malized between themselves at the one-loop level. This 
comes from the fact that the integration of modes in the 
differential slices of energy ±A only produces contribu- 
tions of order ~ dA through particle-particle corrections 
to the above mentioned vertices [Q. The scaling equa- 
tions for Vi n tra and Vumk take the form 

A^fe^ = c n( M - A) (V 2 ntra + V 2 mk ) (2) 



OA 

A = 2c n(fl - A) VintraVumk 



(3) 



with c = 1/^/1 — 4(£'/t) 2 
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(a) (b) 
FIG. 2. BCS vertices that undergo renormalization by par- 
ticle-particle diagrams. The solid and dashed lines stand for 
modes in the neighborhood of the two different saddle points. 
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FIG. 1. Scaling of the chemical potential as a function of 
the high-energy cutoff. The results correspond to the Hub- 
bard coupling U = At. 

The scaling of the chemical potential /Lt(A) towards 
the VHS triggers the different instabilities in the sys- 
tem. These can be traced back to the divergent behavior 
of some of the interaction channels. Among all the kine- 
matics, two different channels develop a singularity for an 
imaginary value of the frequency, namely the BCS chan- 
nel of colliding particles with zero total momentum and 
the channel with momentum transfer Q = (jr, it). Such 
features are related to a phenomenon of condensation and 
the opening of new phases, while the singularities in the 
rest of the channels arise for real values of the frequency 
and they correspond to the appeareance of new excited 
states. 

In what follows, we concentrate on the renormalization 
of the vertices with BCS kinematics, listed in Fig. ||, and 
of those with momentum transfer Q = (7r,7r), listed in 
Fig. |^. In a wilsoninan RG approach of the kind devel- 
oped in Ref. Q| , processes with the latter kinematics are 
disentangled from the divergences in the BCS channel, at 
least at the one-loop level. This holds in particular for a 
saddle-point dispersion relation since, for kinematics in 
which the total momentum of the incoming particles is 
not precisely zero, particle-particle diagrams give irrele- 



If one were to take the bare coupling U of the Hub- 
bard model as the initial point of the RG scaling, the 
setofEqs. (|)-(|) would not provide interesting physics. 
This is so because we would have in such case Vintra = 
V um k — U, and sending the cutoff to zero would simply 
reduce monotonically the interactions in these channels. 
However, the diagonals in this space of couplings mark 
the boundary between the regions of stable and unstable 
scaling. The slightest perturbation with Vintra < V U mk 
will make the couplings to grow large, pointing at the 
appearance of new features in the system. 

kf P+QI k| P+QI 
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(a) (b) 
FIG. 3. Exchange and Umklapp vertices that undergo 
renormalization by particle-hole diagrams. 

An interesting effect comes from the fact that, even 
in the Hubbard model, there are corrections that vanish 
in the low-energy limit A — > (irrelevant operators) but 
may drive to the unstable regime at the early stages of 
the scaling. The Kohn-Luttinger effect, that leads to a 
pai ring instability in the Fermi liquid at very low energies 
]18[ , can be established on the same grounds in the RG 
framework In the case of the Hubbard model, the 
corrections are given by the horizontal iteration of the 
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bubble diagrams in Fig. [|. The momentum flowing to 
the bubbles is not in general close to 0, in the first case, 
nor to Q, in the second. In the wilsonian RG scheme, 
these particle-hole bubbles scale therefore as ~ JdA) 2 , 
and they are considered irrelevant contributions p3. 

In the Hubbard model with t' < 0.276i, the correc 
tions coming from the iteration of diagram (a) in Fig 
have smaller strength than those from diagram (b) 
As long as these are antiscreening diagrams, i.e. they 
add to the bare repulsive interaction, the conditions are 
met to have an unstable scaling for the vertices in Fig. 

with the above constraint on t and t'. The singular 
behavior develops for the combination Vintra — Vumk, in 
which the particle-particle diagrams build up a pole for 
an imaginary value of the frequency. The physical inter- 
pretation is that there is a condensation of Cooper pairs 
at the scale where the vertex function diverges [ H| . The 
symmetry of the order parameter turns out to be d-wave, 
as the combination Vi n t ra — Vumk corresponds to having 
opposite amplitudes in the two saddle points. 
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(a) (b) 
FIG. 4. Particle-hole corrections to the BCS vertices in the 

Hubbard model. 



In our approach, the scaling of the chemical potential 
/x(A) regularizes the divergence of the density of states 
in Eqs. (||)-(||). The result of computing the scale at 
which Vintra — Vumk diverges, for a value of the Hub- 
bard coupling U = it, has been represented in Fig. ||. 
A singularity is only found for values of the bare chemi- 
cal potential /io = A'(Ao) which lead to attraction of the 
Fermi level to the VHS. The curve of the critical scale 
reaches a maximum for a certain value of optimal dop- 
ing, and then it decreases for smaller values of fio as the 
chemical potential is not precisely pinned to the VHS in 
the low-energy limit. 

The instability in the BCS channel has to be matched 
against the strong tendency towards a magnetic instabil- 
ity at wave- vector Q = (n,ir) for t' < 0.276t 0]. In the 
RG framework, this comes from the existence of particle- 
hole contributions which grow large at low energies. They 
are built from diagrams similar to (b) of Fig. 0, but with 
the momentum flowing into the bubble being precisely 
Q. This kinematics corresponds to the vertices Eq and 
Eu that we have defined in Fig. ||. 

It can be seen that, at least at the one- loop level, the 
vertices Eq and Eu are only renormalized by interactions 
with their own kinematics in the wilsonian approach. 
This leads to the scaling equation 



d(E Q + Eu) 
OA 



-c'iEQ+Euf/^H) 



(4) 



In the low- 



with d = log (l + y/l -4(<'/i) 2 ) /(2f/t) 

energy limit, a singularity is reached at a certain value 
of A. The logarithmically divergent particle-hole bubbles 
building up the singularity have an imaginary part equal 
to in/2 times d / (47r 2 i), which means that the vertex gets 
actually a pole for an imaginary value of the frequency. 
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FIG. 5. Plot of the scale of the transition to the pseudogap 
phase (thick line) and of the superconducting instability (thin 
line) as a function of the bare chemical potential and r = t'/t. 

Thus, the divergence in this channel gives rise to 
the condensation of particle-hole pairs with momentum 
Q. This physical interpretation is similar to that of 
the superconducting instability |L9|], but the presence 
in this case of a macroscopic number of particle-hole 
pairs leads to nonvanishing expectation values of the type 
/ d 3 fc(*^ T (k)*Bi(k + Q)> = A, * QCT (k) being the elec- 
tron field operator at saddle point a. This has a dras- 
tic effect in the low-energy spectrum, since the electron 
propagator is corrected by the insertion of the condensate 
as shown in Fig. ^. The result is that a gap of magnitude 
A | opens up in the spectrum of quasiparticles. The ap- 
pearance of this gap takes place in the neighborhood of 
the saddle points, while gapless quasiparticle excitations 
still exist in the rest of the Fermi line not affected by the 
nesting of such singular points. 

The point at which the vertex Eq+Ejj diverges can be 
estimated from Eq. fli]), and it has been represented as 
a function of t 1 in FigTTq. In accordance with the known 
results about the renormalization of t' with doping [fl6[ , 
we have assumed a variation of its effective value follow- 
ing the qualitative trend shown in the figure. The scale 
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of condensation of particle-hole pairs is larger than the 
scale of the superconducting instability up to the point 
of optimal doping marked by the peak of the latter. In 
this regime, the opening of a gap of magnitude |A| is 
the effect that prevails, and the renormalization of the 
interactions is stopped at that energy scale. 

A* A 

k k k I k+Q I k 

^(0) ^(0) „(0) „ 

Lr A Lr A Lr A Lr B Lr A 

FIG. 6. Self-consistent equation for the dressed electron 
propagator G in the particle-hole condensate, in terms of the 
undressed propagators G^ A and G^P at the two inequivalent 
saddle points. 



scale of the superconducting instability, which takes over 
for higher doping levels. The properties that we have dis- 
cussed rely on the existence of an attractive fixed-point 
in the scaling of the chemical potential of the 2D sys- 
tem. They are therefore robust as they do not depend 
on fine-tuning or on the particular details of the model, 
what may explain some of the universal properties of the 
hole-doped copper-oxide superconductors. 

Fruitful discussions with F. Guinea are gratefully ac- 
knowledged. This work has been partly supported by CI- 
CyT (Spain) and CAM (Madrid, Spain) through grants 
PB96/0875 and 07N/0045/98. 



On the other hand, the superconducting instability 
takes over beyond the optimal doping. The chemical po- 
tential departs then from the VHS at low energies, by an 
amount that becomes of the same order as the scale that 
would be predicted for the magnetic instability. In these 
conditions the scaling of Eq + Ejj, which relies on the 
pinning to the VHS, is arrested before the singularity in 
the vertex is reached. This is in contrast to the scaling 
in the BCS channel given by Eqs. (§)-(|j|) , which holds 
irrespective of the pinning to the VHS H ] . 

A closer look at the origin of the singularity in the 
particle-hole sector shows that it arises without a pre- 
ferred direction of the spin. The divergence in the ex- 
change channel may result in a nonvanishing expectation 
value with spin projection S x or S y , depending on the 
choice of the condensate wavefunction. A similar diver- 
gence in the channel with momentum transfer Q and no 
spin-flip shows that J c? 3 fc(^^ T (k)\l/ B |(k + Q)} has also 
an absolute value precisely equal to |A|. Depending on 
the different combinations of phases, the spin of the con- 
densate may point in any direction of the space, what 
provides a check of the rotational invariance of the wilso- 
nian RG scheme. 

When considering the model strictly at zero tempera- 
ture, the ground state of the system is forced to choose 
a definite projection of the spin. The condensation of 
particle-hole pairs leads then to the spontaneous break- 
down of the 15*0(3) invariance. As a consequence, a pair 
of Goldstone bosons arise corresponding to spin waves 
on top of the condensate of particle-hole pairs. In the 
regime below optimal doping where the magnetic insta- 
bility prevails, these are the gapless excitations of the 
spectrum together with the quasiparticles from the re- 
gions not affected by the nesting of the saddle points. 

We have thus clarified the nature of the magnetic in- 
stability that arises when the Fermi level is pinned to the 
VHS of the electron system. The main physical effect is 
the condensation of particle-hole pairs, which results in 
the opening of a gap for the quasiparticles in the neigh- 
borhood of the saddle points. This is the dominant insta- 
bility up to the optimal doping marked by the peak of the 
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